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numerical details of our calculations. 
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1. Introduction 

The b quark mass enters in the theoretical predictions of B mesons decay rates, which provide 
some of the most precise constraints in the unitarity triangle analysis. The Particle Data Group [§] 
quotes a value of rrt^im^) between 4.1 and 4.4 GeV in the MS scheme, as the result of different 
determinations which include lattice computations. 

On the lattice, heavy quarks are treated by using effective theories, in particular HQET when 
dealing with heavy-light systems. The b quark mass can then be predicted using the mass of 
the B meson as experimental input. In this framework, the bare parameters must be determined 
non perturbatively, in order to have a well defined continuum limit. This programme has been 
completed for the lowest order of the HQET (static theory) in the quenched approximation. The 
result from [9], m^m^) = 4.12(11) GeV is in reasonable agreement with [Hj, suggesting that the 
l/m corrections to the b quark mass are small. Here we explicitly compute these corrections 
following the strategy proposed in |JXjj. This clearly represents an important test of HQET, which 
provides indications about its "convergence". 

The method relies on matching QCD and the effective theory in a small volume. The observ- 
ables are then evolved to large volumes where contact with experiments can be established. Such 
an evolution is performed in the effective theory. Different systematics errors have to be controlled 
in the two regimes. For a precise matching in the small volume a very accurate determination of 
the relevant renormalization constants in the full theory is needed. On the other hand, in large 
volume simulations the main problems usually concern the separation of the ground state from the 
excited states.We present details on the numerical evaluation of the various ingredients entering the 
computation of the l/m corrections to the b quark mass. The increase in the (numerical) effort is 
moderate compared to the static case. A related investigation is described in [Qj. 

2. Finite volume 

2.1 Relativistic part 

The observables <l>i 2 are defined as in eqs. 1 (1.2.5), (1.2.6) by (T = L/2) 

4>i(L,M)= <$>\ el (L,M) + <^f it (L) and <J> 2 (L,M) 

where <K el (L,M) = In ( 11^1^] a nd <£f at (L) 

\f\ (0 ',T) J 1 V ; 

The observables <J>j el and <I>2 are computed in (quenched) QCD, with the 0(a)-improved Wilson 
action as in The light quark mass is fixed to be zero. While we use different resolutions 
L\/a = 20,24,32,40, the physical size of the lattice is kept fixed by determining the bare coupling 
such that the SF coupling is 2 g 2 (L[/4) = 1.88 [Q, §]. Using ro = 0.5fm to convert to physical 
units, one finds L\ w 0.36 fm. We take three different values 6,6' = 0, 0.5, 1 for the periodicity 
phases. We want to interpolate our data to the b quark mass, and since according to [J3j] , the value 

'We refer to the equations of [|l|] with a prefix "I". 

2 We have not yet taken the error in g 2 (L\ ) into account. It is expected to be negligible compared to the other errors. 
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Figure 1: Continuum extrapolation of <t>2, for L\M = 12 

of L\M^ Gl is around 12, we simulate three different values of z = L\M : 10.3, 12, 13.2. For each 
of these values <J>2(Li,M) is extrapolated to the continuum via a linear fit in {a/L\) 2 using only 
the lattices of resolution L\/a = 24,32,40. We confirmed that adding the coarsest lattice does not 
change the values of the continuum limit, as one can also see in fig. [I]. We then fit the continuum 
values - see also eq. (1.5.5) - to the following form 



<t> 2 (L u M) = Const + S xLjM, 



(2.3) 



to find the slope 5 = 0.61(5). At this point, we note that a source of error originates from the 
determination [||] of the (u 
percent uncertainty on <I>2. 



determination [[J of the (universal) renormalization factor M^ GI /msF(l/^o)- It amounts to a one 



2.2 Effective theory 

Using the static actions denoted by HYP1 and HYP2 in [||], we compute boundary to boundary 
correlation functions of the QCD Schrodinger functional, but now in the effective theory and for 
two values of L. The smaller one is again L = L\, but with the resolutions L\ja = 6,8,10,12, 
corresponding to j8 = 6.2204, 6.476, 6.635, 6.775 respectively and the larger one is L = L 2 = 2L\ 
with the same values of a. As we did for the relativistic part, for each of these volumes and 
resolutions, the simulations are done with three values of the time extent, namely T = L/2, L/2 — 
a, L/2 + a, and three values of the periodicity angles. Following [[]]], we compute /j stat and /j™ 1 
for these different sets of parameters and from those the functions r s j tat , r'p" and /?^ n as they are 
defined in eqs. (1.2.8 - 1.2.10) - at finite values of the lattice spacing. 

3. Step scaling functions 

Once we have the quantities described in the previous section, we compute the step scaling 
functions a m , a^" , defined as in eqs (1.3.1 - 1.3.3). At finite lattice spacing, we have 

(3.1) 



(u,a/U) = 2U [rf t (2L 1 ,a)-rf t (L 1 ,a)] ii 



T}t«{u,a/L x ) 



Rf n (2L u a) 



Rf n (L u a) 



u=g 2 {U 



Ef(„, fl /L l ) = 2Z., r - n ' 2L '-°»- r " n ' i -°> 
2 V 1 ' R\ m (L u a) 



■ (3-2) 



u=f(U) 
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Figure 2: Continuum extrapolation of the static step scaling function 

Their continuum limits 

o m {u)= lim I m ( M ,a/Li) off (w) = lim T$(u,a/L{) (3.3) 

a/Li^0 ' a/Li^0 

are well defined. Since the static theory is 0(a)-improved, G m (u) is obtained by a linear fit in 
(a/Li) 2 , and we choose to use only the two finest lattices, corresponding to L\/a = 10, 12 . The 
kinetic step scaling functions off (u) are extracted by fitting the data of the three finer lattices 
linearly in a/L. Both for Zf" and for £ m we could have included one coarser lattice without a 
significant change (apart from smaller errors) as illustrated in figs. ^[ 3 

Since we find well compatible results from the different actions after taking the continuum 
limits, we choose a constrained fit in the final analysis. In that case, the error is computed by a 
jackknife procedure, binning the data in such a way that the number of jackknife configurations is 
the same for the different sets of data. We note that a very significant 6,6' dependence is still left 
in the kinetic step scaling functions (1.2.8), nevertheless the final physical results will turn out to 
be independent of these kinematical variables as they should (apart from small l/m 2 terms). In our 
figures, we made the choice (d,d f ) = (0,0.5). 

4. Large volume 

The quantities computed in "infinite" volume are the static and the kinetic energy (1.4.2), where 
now the light quark mass is fixed to m stom g e . These energies occur in eq. (1.4.1), through the terms 
[£ stat - r S ! tat (L 2 )] and [o^ n (£ kin - rf n (L 2 ))} . The infinite and the finite volume part of the static 
and the kinetic energy are of course computed at the same value of the lattice spacing (in order 
to cancel the divergences, as we did for the step scaling function). Our approximation to infinite 
volume is Leo = AL\ ~ 1.5 fm (again using ro = 0.5fm), and we use two resolutions L^/a = 16,24, 
corresponding to /3 = 6.0219, 6.2885 respectively. We obtain a nice signal for E stBt due to the 
HYP 1/2 actions ||]. The errors in the effective mass remain below a level of few MeV out to 2 fm 

3 Since the definition of £,„ and E^™ involves simulations on a lattice of time extent T = L/2 — a, a space extent of 
L\/a = 6 is too small to obtain trustworthy numbers. For this reason in the plots of these step scaling functions, we do 
not show the data coming from this lattice. 
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Figure 3: Continuum extrapolation of the kinetic step scaling functions 

time separation. Thus the significant excited state contamination can be removed with confidence 
and £ stat -rf at (L 2 ) = 612(2) MeV and 597(2) MeV are found at j8 = 6.0219, 6.2885 respectively 
for the action HYP2. From this we quote at present £" stat — Pj tat (L2) = 586(1 l)MeV as a continuum 
number. 

However, [o\: m (E km — Tj m (L2))] is very small and the errors of the effective (xo-dependent) 
2? kin g row ra ther rapidly with xq. Thus, although we will see that this error is not dominant in the 
end, [(OianiE^ — r^ n (L 2 ))] has a large relative error. 



5. The RGI b quark mass 



5.1 Static part 



The value of the RGI quark mass in the static approximation, defined by eq. (1.5.4), is obtained 
by a linear interpolation of (1.5.2), where we used m e w enmental _ 5404 MeV as appropriate for a 
strange quark. This interpolation is illustrated in fig || We obtain (the conversion to physical units 
is done with rg = 0.5 fm) 



roMt 



17.18(25) (HYP2) 
17.15(25) (HYP1) 



J 6771(99) MeV (HYP2) 
[ 6757(99) MeV (HYP1) 



(5.1) 



5.2 l/m corrections 

The next to leading order correction of the b quark mass can be separated in two parts, cor- 
responding to wig , computed entirely from the small volume, and mg lb \ given by eq. (1.5.3). 
Following eq. (1.5.6) we have 



M 



(la) 



Ly))<t>y(L U Mr) 



i lb) = -^[^(£ kin -r^(L 2 ))^ n ], 



(5.2) 
(5.3) 
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Figure 4: Interpolation and solution of eq. (1.5.2) 



and we find for the action HYP2 
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-0.08(4) 


-0.14(11) 
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-0.08(4) 
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-0.08(4) 
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14 



while for HYP1 the numbers are compatible but have larger uncertainties. It is now clear that the 
results are 6, ^'-independent. In physical units, we find the scale and scheme independent numbers 



(la) 



-30(15) MeV Af, 



(lb) 



-56(43) MeV 



(5-4) 



in the quenched approximation. 



6. Conclusions 

Our numerical results show that indeed the l/m corrections can be computed with a precision 
which is - at present - better than the one of the leading (static) result: the absolute errors of leading 
and next-to-leading terms are in a ratio of about two to one. This is a first demonstration of the 
practicability of the strategy of ^ beyond the leading order. Still, as mentioned in the previous 
section, the errors of correlation functions with the kinetic insertion grow rather rapidly in time and 
when decreasing the lattice spacing. The former behavior is at present not understood. 

The size of the l/m corrections M^ la ' and M^ lb ' is small and confirms the validity of the naive 
counting in terms of powers of A/mt,. Further numerical confirmation of this is provided by the 
close agreement of our result with the one - also in quenched QCD - of [|J], r^M^ = 17.08(41) where 
the effective theory was not used. This also means that l/m 2 terms can be neglected altogether. 

On the other hand, our static result differs (statistically) significantly from the static one of [0], 
rnMb — 16. 12(25) (15), where the second error is in common between that computation and the 
present one. In ^ the matching was performed using an other observable, and at L = Lq = L\/2. 
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By a comparison to our present NLO result, we see that this leads to a somewhat larger (~ 5%) 
l/m correction. This is not at all unexpected, cf. item 3 in section 6 of [1]. In fact, an explicit 
computation of this correction would be another interesting confirmation of the applicability of the 
effective theory 

Let us finally translate our numbers to the MS scheme. With ro = 0.5 fm and = 238MeV [0] 
we find for the b quark mass at its own scale 

m h = mf lt + m£ ) , mj; tat (mb) =4.350(64) GeV, m£\m h ) = -0.049(29) GeV. (6.1) 

Despite the employed quenched approximation, the total result, m\,(m\,) = 4.30(7) GeV, is in 
good agreement with the range quoted in the particle data book and not far from the precise value 
of m b (mb) = 4. 19(5) GeV derived in [[J from the e + e~ — > b total cross section and high order 
perturbation theory. Our result is also compatible with the early computation (quenched and static) 
of Jl0| ] m\,(m\,) = 4.41(5)(10)GeV, where the matching was done perturbatively at the NLO. This 
last result was updated to 4.30(5)(5) GeV in [[H]] with the help of the NNLO result of [@|. 



It seems that on the one hand it is the right time to apply these methods both to quantities 
where larger l/m corrections are expected, such as Fg, and to the mass of the b quark in full QCD. 
On the other hand, research should continue to improve the statistical errors in particular of the 
l/m corrections in large volume. A promising route is to follow [|l3|, [i4| ]. 

Acknowledgement. We thank NIC for allocating computer time on the APEmille computers at 
DESY Zeuthen to this project and the APE group for its help. 
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